Vector meson photoproduction studied in its radiative decay channel 
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We provide an analysis of vector meson photoproduction in the channel of the vector meson 
decaying into a pseudoscalar meson plus a photon, i.e. V —> P<y. It is shown that non-trivial 
kinematic correlations arise from the measurement of the P7 angular distributions in the overall 
cm. system in comparison with those in the vector-meson-rest frame. In terms of the vector 
meson density matrix elements, the implication of such kinematic correlations in the measurement 
of polarization observables is discussed. For the u) meson production, due to its relatively large 
branching ratios for lu — » 7r°7, additional events from this channel may enrich the information about 
the reaction mechanism and improve the statistics of the recent measurement of polarized beam 
asymmetries by the GRAAL Collaboration. For <f> — » 777, p — » 717, and K* — > K'y, we expect 
that additional information about the spin structure of the vector meson production vertex can be 
derived. 

PACS numbers: 24.70.+S, 25.20.Lj, 13.88. +e, 13.20.-v 

I. INTRODUCTION 

Direct experimental evidence for baryon resonance couplings into vector meson has been searched for in vector 
meson photoproduction reaction in recent years Q, 0, 0, Q ■ One of the motivations behind this effort is to find 
"missing resonances" , which are predicted by the nonrelativistic constituent quark model (NRCQM) but not found in 
ttN scatterings 0, @ • The study of vector meson photoproduction at large angles near threshold is particularly useful 
due to the relatively small contributions from the background processes, such as £-channel natural and unnatural 
parity exchanges. 

Various experimental projects have been and are being carried out at ELSA, JLab, ESRF, and Spring-8. In addition 
to the cross sections, polarization observables, which are more sensitive to resonance excitations, will also be measured 
using newly developed experimental techniques 0, @ ■ This is of great importance for the purpose of disentangling 
the s-channel resonance excitations in the reaction. Since a large number of degrees of freedom will be involved and 
due to the lack of dynamic information in such a non-perturbative region, one, in principle, needs a complete set 
of measurements of all independent spin polarization observables to obtain sufficient information about resonance 
excitations and their couplings to the meson and nucleon 0. 

The polarization observables and density matrix elements are essentially equivalent languages that connect the 
theoretical phenomenologies with the experimental observables. The density matrix elements, as the interference 
between the independent transition amplitudes, can be measured in the final state vector meson decay distributions. 
They can be directly compared with theoretical calculations, and hence give access to dynamical information about 
the transition mechanism. For the vector meson photoproduction reaction, A number of analyses have been carried 
out in the literature. Density matrix elements for polarized photon beams were discussed by Schilling, Seyboth 
and Wolf 10J. More recently, Pichowsky, Savk h and Tabakin studied various relations between polarization 
observables and helicity amplitudes. In Ref. IT^. Kloet et al. investigated inequality constraints on the density 
matrix elements defined for the produced vector meson with unpolarized or linearly polarized photon beams. They 
also explored relations between vector meson decay distributions in the vector-meson rest frame and the overall 7-JV 
frame, which can be extended to the analysis of other decay channels. 

Another aspect relevant to the derivation of the density matrix elements is the dynamics for the produced vector 
meson coupling to the detected particles. For instance, u> meson is usually detected in u> — > 7r + 7r~7r° and p° in 
p° — > 7r + 7r~. The angular distributions for vector meson decays into spinless particles were studied in Ref. [T(J- 
However, their decays into photon and pseudoscalar meson, V — > P7, have not been discussed. In this work, we will 
develop the formalism for vector meson decay into photon and pseudoscalar meson, e.g. u> — * tt°7, 4> — > fff, p — ► 7r7, 
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and K* — > ^7. We will show that additional information about the vector meson production mechanism can be 
obtained from measurements of those final states. This analysis has an advantage in u meson photoproduction, of 
which the branching ratio of uj — > 7r°7 is about 8.5% 0] . As an isoscalar meson, only nucleon resonances can contribute 
in its s- and w-channel production. This reduces significantly the number of excited states in the corresponding 
reaction channels. In this sense, the additional information from uj — ► 7r°7 should be useful for constraining the 
nucleon resonances in — ► ujN. Nevertheless, since the uj decays are dominantly via uj — ► 7r + 7r~7r°, and ir + ir~ 
could form a configuration of J p = 1" as a photon, it is interesting to compare these two distributions and gain some 
insights into the sensitivities of the vector meson production mechanism to these two decay channels in a polarization 
measurement. This possibility is enhanced in the comparison between V — > PP and V — > P7. For 4> — > 777, the 
branching ratio is still sizeable. In comparison with the dominant decay channel of <p — > KK, the decay distribution 
of — » 777 contains different spin structure information; hence should be more sensitive to the production mechanism. 
Similar feature applies to the p — ► 7r7 and K* — ► K"/ in p and K* meson photoproduction. 

In this paper, the analysis of the decay channel V —* P7, will be presented in Sect. II. In Sect. Ill, we discuss 
the kinematic correlations arising from this decay channel and analyze their influence on the measurements of spin 
observables in the the overall cm. system. Comparison with the results derived in the vector-meson-rest frame 
will then be made. A summary will be given in Sec. IV. To make it convenient for readers, analyses of the spin 
observables in terms of the density matrix elements of the polarized particles are included in the Appendix. 



II. DENSITY MATRIX ELEMENTS FOR V — > P7 

In the overall cm. system, the invariant amplitude is defined as 

T KXf . XiXi = (q,A, ) ;P / ,A / |r|k ) A 7 ;Pi,A i ) -» (A„A/|T|A 7 Aj), (1) 

where k and q are momenta of the initial photon and final state vector meson, respectively; Momenta Pi = — k and 
Pf = — q are for the initial and final state nucleons, respectively; A 7 (= ±1), A u (— 0,±1), Xi (= ±1/2) and A/ 
(= ±1/2) are helicities of the photon, vector meson and initial and final state nucleons, respectively. 

The amplitudes for the two photon polarizations A 7 = ±1 are not independent of each other. They are connected 
by the Jacob- Wick parity relation: 

(A„A / |T|A 7 A l ) = (_l)(^-A/)-(A,-A i)( _ At; _ A/ | T | _ ^ _ A .). (2) 

We first consider the vector meson decay in the overall 7- TV cm. system. The decay of vector meson (J PC = 1 — ) 
into a pseudoscalar meson (J p = _ ) and photon (J PC = 1 ), i.e. V — > P7, is described by the effective Lagrangian 

£ V p 7 = ^L €a0 sd a A^d^V 5 P , (3) 

My 

where gy is the coupling constant and e Q( a 7 5 is the Levi-Civita tensor; V s and P here denote the vector and pseu- 
doscalar meson fields, and A@ is the photon field; This effective Lagrangian is the only Lorentz-covariant one for the 
VP7 coupling at leading order. 

We assume that the fields couple as elementary ones. Therefore, the momentum-dependence of gy can be ne- 

1/2 

glected, i.e. gy = gy(0). This assumption depends on the average radius (r^) of the vector meson, for which we 
phenomenologically refer to the average disctance between the quark and antiquark within the meson. As a simple 
estimate for the ground state mesons of which the spatial wavefunctions are spherical, the particle size effects appear 
as the quadratic term in the form factor 1 — (r^)/6, where k q ~ 380 MeV are the typical constituent quark momenta. 

1 /2 

For a typical size of (r|) ~ 0.5 fm for the mesons, a correction of about 15% seems to be needed. However, this will 
not affect the general feature of the subsequent analyses due to two reasons: i) As shown by Becchi and Morpurgo [T^| . 
in the limit of k q — 0, i.e. gv(k q ) — gv(0), the calculations for both limits of My = Mp{exp) and Mp — My(exp) 
give the same results which are consistent with the experimental data. This suggests that the momentum dependence 
of gy is not significant, ii) Even though gy(k q ) ^ ffv(0), it will only lead to some percentages of descrepancies in 
the event countings between the calculation of Eq. J3Jl and experimental statistics, while the decay distributions will 
not change. In this sense, we can neglect the the momentum dependence of the gy at this moment, but focus on the 
study of kinematic correlations arising from the measurement of the P7 angular distributions. 

It is worth noting that the vector meson size effects may become non-negligible at the production vertices for 
7./V — > VN, and careful considerations of the vector meson coupling form factors are crucial. However, this is 
independent of our motivation here. As mentioned in the Introduction, we are interested in possible experimental 
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measurements of the angular distributions of final state particles, through which dynamical information about the 
production mechanism (contained in the density matrix elements) can be extracted. 

As illustrated in Fig.^ we define the Zf axis along the three moment of the vector meson, the y/ axis is normal to 
the production plane defined by k x q, and hence the Xf direction is determined by yf x Zf. Since the vector meson 
three momentum |q| is determined by the cm. energy W , the three momentum of the final state photon in the overall 
cm. system can be expressed as p = q/2 + r, where r is the three momentum of the photon in the vector meson rest 
frame. The z / axis also defines the z direction of the decay co-ordinates, while the polar and azimuthal angles 9 and 
4> are defined by the momentum difference between the final-state photon and pseudoscalar meson, p — p s = 2r. On 
the other hand, for the description of vector meson decay, it is convenient to select the decay angles, 9 C and (j) c , which 
are the polar and azimuthal angles of the flight direction of the photon in the overall cm. system with respect to Zf. 
As shown in Fig. |3 we have, 



posmfc = |r|sinc/, 
Pocos9 c = |r|cos# 



where _Po = |p| is the energy of the final-state photon. 

The decay amplitude of the vector meson with transverse polarization can thus be expressed as: 



(4) 



(A 7 ;# c ,0 c |M|A„ =±1) = C\ 



3 p (g°A 7 + |q|A„) 
M v 



8tt 

and the amplitude with longitudinal polarization is: 

(A 7 ;# c ,0 c |M|A.„ =0) = C\ 

= C\ 



c HT p (g°A 7 + |q|A,,) , ^K D i* 



M 



— |p|A 7 e 7 q 
'^IP|ArDo%( 



(5) 



(6) 



where A 7 = ±1 and At, = 0, ±1 are the helicities of the photon and vector meson in V — > P7, respectively. The 
coefficient C = egy is the coupling constant. 

The Wigner rotation functions follow the convention of Rose [l5| : 



D I MN {a,P, 1 )=e-* Ma + N ^d I MN {(3), 



(7) 



where a, (3, and 7 are Euler angles for the rotations of a vector. 

Similar to the case of vector meson decaying into spinless particles [Tol | , we can express the angular distribution of 
V — * P7 in terms of the vector meson density matrix elements p\ v \> : 



dN 



dcos 



W(cos c ,(f> c ,p) 

Wtt(cos0 c ,(/> c ,p) + W LL {cos9 c ,(j) c ,p) + W T l(cos0 c ,<I>c,p), 



(8) 



where 9 C and ip c are functions of 9 and 0; and the latter ones are the variables of the differential distributions 
on the left-hand side; the subscription TT , LL and TL denote the interfering distribution functions between the 
transverse-transverse, longitudinal- longitudinal and transverse- longitudinal vector meson polarizations. 
They have the following expressions in the overall cm. system: 



Wtt{cos0 c ,4>c,p) = — 



A„,A 

C^3_ 

(To 87T 



E 

=±1:A,,A1 



(A 7 ;^0|M|A„) /9A „A^ A A,(A;|Mt|A; 



E 

A„,A'„=±1:A T 



p (?oA 7 + |q|A„ 



-(-1) A ^-A„A 7 ( 



-<l>c)p\ v \> v (V) 



C 2 3 vl 



p (g A; +]q|A4) 
' Mv 



cr 8tt 2 



{ [Tl{9 c ) + J%(0 C )] (pu + + sin 2 9 C [e~ 2 ^ + e 2 ^"p-ii] } 



(9) 
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and 



C 2 3 

W LL (cos6 c ,<f> c ,p) = — — V |p| 2 A 7 A'r^ Uc,O c ,-ct) c )S A A'Poo(V)Dl A , (0 C , 6> c , -0 C ) 
an air * — ' ' 7 ~ t 



fin 87T 

A„=A'„=0;A 7 ,A^ 



— i-|p|2 sin 2 0cpoO; (10) 
(Tq 07T 



and 

CP3_ po|p|A . 

cr 8?T ^ My 



W Ti (c08fl c ,^,p) = — — ^^{(9oA 7 + |q|AO(-l) A "I?^A 7 (^^ C ,-0c)pA„O^A T (0c^c,-0 c ) 

ff0 87r A.,Al^0;A 7 MV 1 

+ (g o A 7 + |q|A;)(-l) V "^(0 c ,0 c ,--^) m , J Di A , A7 (^^ c ,-0 c )} 
( ~ " /jr,H "^-(^i(^)-^2(ec))[ e -^(pio-po-i)+e^(poi-P-io)], (11) 



where 



Fi(6 D ) = r-L(®> + |q|)(l-cos0 c ), 

My 

^2(9 C ) = ^-((Zo-|q|)(l+co S c ) (12) 

are functions of C1 and hence introduce kinematic correlations to the coordinate transformation. The 5 function 
#A 7 A' implies the sum over the final-state unpolarized photons. The factor ero = C 2 |r| 2 is the normalization factor 
for V — ► Pj, and is proportional to the decay width for V — > P7 in the vector meson rest frame. Note that 
VF(cos # c , 4>ci p) is coordinate-dependent due to the Lorentz transformation from the vector meson rest frame to the 
overall cm. system. 

The density matrix elements of the vector meson are independent of the co-ordinate frame selection for the vector 
meson decays. They are related to the initial photon polarizations in the overall cm. system via: 



'A V ^ = 7J E ^„j /! i 1 A i /'A 1 j;(7)^A / i;j, , (13) 

A/A 7 A;A^ 



PKK, ■ N 



where N = \ ^2 X ^ \ x I^a„A/,a t A; | 2 is the normalization factor and is double the unpolarized cross section ignoring 
the phase space factor. 

The polarized-photon density matrix element ^0(7) is defined as |ldj |: 

p(H) = \{I 1 + CT-V 1 l (14) 

where er is the Pauli matrix for the photon's two independent polarizations; P 7 determines both the degree of 
polarization (via its magnitude P 7 ) and the polarization direction. For linearly polarized photons, with $ denoting 
the angle between the polarization vector of the photon (cos^, sin<I>, 0) and the production plane (x%-Zi plane), one 
has P 7 ^~ Pf 

(— cos 2$, — sin 2$, 0). For the circularly polarized photons, the polarization vector is along the Zi axis 
and hence P 7 = P 7 (0, 0, A 7 ) with A 7 = ±1. 

Substituting Eq. (|14|l into H13[l , we obtain the familiar form of the vector meson density matrix elements |l(1 | : 

Pa„AJ, — '^T7 E -^t^/^T^i-^A^Aj.A^AiJ 



27V 

A-y A f A 



■2 



^ \ ^ rp rji* 

2^ 2^ Ja - a / --Mi J A'„A / ,A 7 A i > 
A-,, A f Xi 

Px v X' v = E ^7^„A / ,-A 7 A i r\^A / ,A- ( Ai' 



I 

A-y A f A 



1 

2N 



^ \ 71 A l ,A / ,A 7 A i T , A;A / ,A 7 A i ■ (15) 
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The decomposition of the vector meson decay distribution in terms of the initial photon polarizations thus allows 
us to express Eq. © as 

3 

W(cos9 cl (j) cl p) (cos 9 C1 ^ p°) + Y,P"W a (cos 9 c ^ c , P a ) , (16) 

a=l 

where W° denotes the distribution with the unpolarized photons, and W 1,2 denote those with the linearly polarized 
photons, and W 3 with circularly polarized photons. 

With the Jacob- Wick parity relation [Eq. J2J] and the requirement that the density matrix elements must be 
Hermitian: 

pZk=pk^> ( 17 ) 

the decomposed distributions can be obtained: 

W°(cos0 c ^ c ,p°) = A^i{sin 2 ^ + ^[^i(^)+^ 2 (^)]Pii 

+ sm 2 e c cos2(j) c p\_ l + V2[T 1 (e c )- T 2 (O c )]sme c cos(l) c Rep w ^ , (18) 

87T (Jo ^ I 

+ sin 2 6» c cos2(/ )cj o}_ 1 + V2[J"i(0 c ) - J" 2 (^ c )]sin6» c cos( ? !) c Rep 10 | , (19) 
W 2 (cos 9 c ,<j> c ,p 2 ) = ^ { sin 2 9 C sin 20 c Imp?_ 1 - v^[^i(0 c ) - -F 2 (0 C )] sin C sin <^ c Imp 2 [ (20) 

07T (Jo I. J 

and 

W 3 (cos0 c ,(/> c ,p 2 ) = J sin 2 c sin20 c Im /9 ?_ 1 - v^[^i(0 c ) - ^ 2 (0 C )] sin0 c smcf> c lmp\ \ . (21) 

07T (Jo ^ J 

Recalling again that 9 C is a function of 0, strong kinematic correlations have been embedded in the above expressions. 
Naturally, one would expect that in the limit of |q| — * 0, and hence q — > My, the azimuthal angles (9 c ,4>c) Wlu be 
identical to (9, </>), and C 2 p 2 /<jq = Pq/\t\ 2 = 1. These expressions then reduce to the ones derived in the vector- meson- 
rest frame |l0|: 

W°(cos9,(t>,p Q ) = {sin 2 0p° o + (1 + cos 2 0)p° n + sin 2 9 cos20p?_ 1 

+v / 2sin26»cos0Rep? o } » ( 22 ) 

3 

W 1 (cos9,(/),p 1 ) = — {sin 2 ^ + (1 + cos 2 6») / 9} 1 +sin 2 6»cos20p}_ 1 
8tt 

+v / 2sin26»cos0Rep lo } , (23) 
W 2 (cos 9, cf),p 2 ) = ^- (sin 2 6»sin20Imp 2 _ 1 + a/2 sin 29 sin 01mp 2 o } (24) 

and 

W^cosfl^p 3 ) = — (sin 2 6»sin2(?!)Imp5_ 1 + V2 sin 26> sin (j)lmp 3 1Q } . (25) 

As shown in Eq. (J5J, the non- vanishing three momentum of the vector meson in the overall cm. system introduces 
an additional term for the transverse vector meson decays. As a consequence, kinematic factors T\ i(Q^), which are 
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functions of the decay angle 9 Cl appear in the TT and TL decay distributions. Nevertheless, co-ordinate transfor- 
mations also exist between 9 C and 9. In comparison with Eqs. I|22l) - (|25[l for vector meson decay in its rest frame, 
kinematic factors in Eqs. H18(l - (|21|l cannot be cleanly factorized out in the measurement of V — > P7 in the overall cm. 
system. This is different from the measurement of vector meson decay into spinless mesons e.g. p° — > n + n~ or 
4> — ► K + K~ , where a kinematic factor can be taken out of the distribution functions by selecting the polar angle 9 as 
the angle between the momentum difference of the final-state pions and the vector meson momentum in the overall 
cm. system. Such a correlation also reflects the dynamical difference between V — > P7 and V — > PP from which we 
expect to learn more about the vector meson production mechanism via the measurement of the decay distributions. 

In the next section, we will discuss some features arising from the kinematic correlations and their impact on the 
polarized beam asymmetry measurement in different frames. 



III. KINEMATIC CORRELATIONS IN POLARIZED BEAM ASYMMETRY 

In the production plane, the measurement of the linearly polarized beam asymmetry is defined as the cross section 
difference between polarizing the photons along the Xi = 0°) and ?/i-axis (i> = 90°), which correspond to P 7 along 
respectively. The cross sections for these two polarizations thus can be expressed as: 

W x ($ = 90°,p) = / / d£W o (cos0 c ,0 c y)-P 7 / / dfW^cosflc^p 1 ) 

= ^V) + ^V) , (26) 

and 

W\\($ = 0°,p) = W°{p°) - P 1 W 1 {p x ) . (27) 
The linearly polarized photon asymmetry is thus defined as: 



g W^ = 90°,p)-W\i^^0°,p) = W\p\ 
W±($ = 9Q o ,p) + W\\(<f> = 0°,p) 7 IT (/9°) ' 



where P 7 is determined by the experimental setup, and W 1 (p 1 ) and W°(p°) are the integrals over (9,4>): 
W\p 1 ) = [ dUW 1 (cos 9 c ,<t> c ,p 1 ) 

dn{[(p ■ qf - M*(p ■ e°)> 00 + [(p • qf + MUP ■ 0>n} , (29) 



and 



3 1 

8^ |r| 2 M£ 



w°(p°) = I dnw° (cos e c ,<t> c ,p°) 

3 1 



8tt |r| 2 M 2 



dn {[(p ■ qf - M 2 v {p ■ e° v f}p° 00 + [(p ■ qf + ■ 0>n} , (30) 



where the relation p 2 [F 2 (9 c ) + T 2 {9 c )]/2 = [{p ■ qf + M 2 (p ■ e°) 2 ]/M 2 has been applied, and e° = (\q\,q q)/M v is 
the longitudinal polarization vector of the vector meson. 

Note that the density matrix elements p a (a = 0, 1, 2, 3) are independent of 9 and cf>, we define the following two 
integrals: 

W a . Ij^f^p-qf-MUp.^f] 

Wb = hwip} J dn[iP ■ 9)2 + Ml{v ■ e ° )2] - (31) 

Hence, the polarized beam asymmetry can be expressed as 

y = POO + (Wb/WJfa 

p° 00 + (W b /W a )p n ' [6 1 
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and the ratio of the integrals gives: 



W b _ 1 |q| 2 L 3 M 



^ = 2+^(1 + ^] , 03) 



where in the integration, the relation 2qopo = + p\ — M P — (r — q/2) 2 has been used. In the vector-meson-rest 
frame, |q| — ► 0, we have Wb/W a = 2, and the polarized beam asymmetry reduces to 

y _ POO + 2 Pll (oa\ 

which is the familiar result derived in the vector meson decays into pseudoscalar mesons |17| . 

Due to the non-trivial kinematic correlations arising from the decay distributions, the polarized photon beam 
asymmetry is also accompanied with a kinematic correlation factor when transforms from the vector-meson-rest 
frame to the overall cm. frame. Although one, in principle, can select the vector-meson-rest frame as the working 
frame for deriving the polarization asymmetry, in reality, it will depend on how well the vector meson kinematic is 
reconstructed. For narrow states, such as u) and <p, the transformation of the working frame from the overall cm. 
one to the vector-meson-rest one should have relatively small ambiguities. However, for broad states, such as p, the 
sizeable width will lead to uncertainties of determing the three moment q. The kinematic correlations hence will 
produce significant effects in the measurement of the polarization observables. 

This kind of situation makes Eq. i|34|) useful for evaluating the kinematic correlation effects. Firstly, note that the 
second term in Eq. I|34|) explicitly depends on the mass and momentum of the pseudoscalar meson in the vector- 
meson-rest frame. It shows that the kinematic correlations will become more significant if the vector meson decays 
into a heavier pseudoscalar meson instead of a lighter one due to the increasing ratio Mp/\r\ 2 . For instance, for the 
final state decays of co — ► 777, M^/\r v \ 2 = 0.13, and for to — > 7r°7, Af 2 /|r ?) | 2 = 7.56. This ratio is essentially a constant 
for a fixed decay channel. Therefore, Eq. (|34l) suggests that above the vector meson production channel and at a 
fixed production energy kinematic correlation effects should show up increasingly significant for heavier pseudoscalar 
decay channels. 

Secondly, for a fixed pseudoscalar decay channel, Eq. (|34[1 shows that the kinematic correlation should also become 
increasingly important with the increasing reaction energies due to the term proportional to |q| 2 /My. In another 
word, only when the vector meson is produced near threshold, i.e. |q| 2 /My << 1, can the kinematic correlation effects 
be neglected as a reasonable approximation. Otherwise, uncertainties from such correlations should be cautioned. 



IV. SUMMARY 



We have carried out an analysis of vector meson decay into a pseudoscalar meson plus a photon in photoproduction 
reactions, e.g. uj — > 7r°7, tfi — > 777, p — > 717, and K* — » Kj, etc. Comparing with the previous measurements of 
vector meson decays into pseudoscalars, this channel has advantages of providing additional dynamical information 
about the vector meson production mechanism due to the different spin structures carried by the differential decay 
distributions. 

It was found that kinematic correlations would become important in this decay channel in the overall cm. system. 
An explicit relation for the correlations between the vector-meson-rest frame and the overall cm. system was derived 
in Eq.(33), which highlighted the kinematic sensitivities of the polarization observables studied in different reactions 
and different pseudoscalar meson decay channels. It is thus useful for providing guidance for experimental investigation 
of vector meson photoproduction for the purpose of studying nucleon resonance excitations in polarization reactions. 

Although the decay channel V — * -P7 generally has small branching ratios for most vector mesons, it is relatively 
large for the omega meson with & w _^o 7 = (8.5 ± 0.5)% [3. Therefore, for to meson photoproduction, this channel 
will not only increase the experimental statistics, but also provide an independent measurement of polarized beam 
asymmetry which can be then compared with the one measured in u> — > 7r°7r + 7r~. For <f> and p°, their dominant 
decays are into two pseudoscalars. Hence, additional spin structure information can be expected in <j) — > 777 and 
p° — > 7r°7. In particular, br^ rn = (1.295 ± 0.025)% is still sizeable, and makes this channel an important source 
for the </> meson photoproduction mechanism near threshold. Experimental facilities at ESRF (GRAAL), SPring-8 
(LEPS) and ELSA (Crystal Ball) with linearly polarized photon beams and charge-neutral particle detectors should 
have advantages for addressing this issue. 

We also derived the single polarization observables for vector meson photoproduction in terms of the density matrix 
elements for the polarized particles. Those elements can be directly related to the experimental measurement of the 
corresponding angular distributions of the vector meson decays into either spinless mesons or a spinless meson plus 
a photon. For linearly polarized photon beams the angular distribution of the photoproduced-vector-meson decay 
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into spinless particles has been developed in the literature [10| . The polarization observables in terms of the bilinear 
helicity product of the helicity amplitudes have also been discussed widely in the literature 0. In the Appendix, we 
presented the vector meson decay distribution functions in terms of those measurable density matrix elements defined 
for those polarized particles. This should be useful for future experimental analyses at GRAAL, JLab, and SPring-8, 
with polarized targets, recoil polarization or beam-target double polarizations. 
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Appendix: Density matrix elements for polarization observables 

In Ref. |10|, the density matrix elements for polarized photon beams were derived. We shall adopt the same 
convention and derive density matrix elements for other single polarization (i.e. polarized target, recoil polarization, 
and vector meson polarization), and double polarization observables. First of all, we will outline several basic aspects 
of density matrix elements with polarized photon beamSjWhich will be useful for further discussions. 

For the convenience of comparing with other analyses |9j , we also express the invariant amplitudes as the following 
12 independent helicity amplitudes: 





— (A„, Xf — 


+1/2|T|A 7 


= 1,A; = 


-1/2) 




= (X v ,Xf = 


+1/2|T|A 7 


= l,Ai = 


+1/2) 




= \X v ,Xf = 


-1/2|T|A 7 


= l,Ai = 


-1/2) 




= \X v ,Xf = 


-1/2|T|A 7 


= l,Xi = 


+1/2) 



A. Convention and kinematics 



The general form of the angular distribution for the vector meson decay into spinless particles (e.g., u> — * ir + ir ir°, 
p° -> 7T+7T-, cf>^ K+K~) is 



dN 



dcos 



ly(cosM) = ]T(M|Af|A> A „ A; (VO(A;|Aft|0,0) , (36) 



where M is the vector meson decay amplitude and can be factorized as 




{6, 4>\M\X V ) = C\ —Di*(^ 0, -0) ; (37) 



The decay angles, 6 and <f>, are defined as the polar and azimuthal angles of the flight direction of one of the decay 
particles in the vector meson rest frame in the case of two-body decay, e.g. p° — » tt + it~ , <j> — > K + K~ . For three-body 
decay of the vector meson, e.g. ui — > 7r°7r + 7r~, 8 and 4> denote the polar and azimuthal angles of the normal direction 
of the decay plane. The constant C is independent of A^ and determined by the vector meson decay width; the Wigner 
rotation functions D follow the convention of Ref. as defined by Eq. Q . 

Consequently, the angular distribution of Eq. Q3(j|) can be expressed in terms of the vector meson density matrices 

VF(cosM) = |- D xlo(&0,-<t>)pKK( V ) D l>A < t ) > 6 >- < f > ') ' ( 38 ) 
x v \' v 

where p(V) is Hermitian, p\^\' v {^) — Pyx (^)- Meanwhile, since W (cos 6, (ft) is a linear function of p(V), it can be 
decomposed into a linear combination in terms of the polarization status of the particles. 

The above equation is general for the vector meson decays into spinless particles. The vector meson density matrix 
element p(V) can be related to the polarization status of the initial and final state particles in their spin space via 



9 



the transition amplitudes T\ v x f ,\ J \ i - Therefore, for different polarization reactions with the vector meson decays 
into spinless particles, p(V) is the source containing information about the polarization observables. The angular 
distribution, as a function of p(V), hence provides access to the physics reflected by p(V). In experiment, p(V) is a 
quantity that can be derived from the data for vector meson decay distributions, while in theory, it can be calculated 
with dynamical models [H [H il0|2p. 

Since p(V) is related to the polarization of the initial and final state particles in their spin space, this is equivalent to 
saying that for different polarizations of the initial and final state particles there are different density matrix elements 
can be measured. In this sense, we are also interested in the number of minimum measurements which can provide 
the maximum amount of information on the transition mechanisms [9j. 

As follows, we will first re-derive the polarized beam asymmetry in terms of the density matrix elements following 
Ref. |lOj| . and then derive other polarization asymmetries in terms of the corresponding density matrix elements. 



B. Polarized beam asymmetry 

Substituting Eq. (|14|l into (|38[1 , we can easily reproduce the results of Ref. [10( and obtain the angular distribution: 
W (cos 0,0,$) = W°(cos9,(p,p XvK ) - P 7 cos2$iy 1 (cos6»,</>,p^ A , ) 
-P 7 sin 2^W 2 (cos 0, </>, p 2 XvK ) 

+A 7 P 7 ^ 3 (cos0,0,p 3 i , A ,) . (39) 

where 

W°(cos9,^p Q ) = A[Isin 2 + i(3cos 2 0-l)p° o 

— v^Rep^o sm 2# cos cp — p\_ x sin 2 cos 2(f>] , 
W 1 (cos9,<j),p 1 ) = ^-[pl lS m 2 e + p 1 00 cos 2 9 

— V2Rep\ sin 20 cos cj> — p\_ x sin 2 6 cos 2(f)] , 
3 

W 2 (cos0, (f>, p 2 ) = — [V2Imp 2 o sin20sin0 + Im / 9 2 _ 1 sm 2 0sin20] , 

47T 

3 

W 3 (cos0,0,p 3 ) = — [V2Rep 3 x sm20sin0 + Imp 3 ^ 1 sin 2 0sin20] . (40) 

47T 

For the polarized photon beams, the density matrix element p° corresponds to the unpolarized photon measurement, 
p 1 and p 2 correspond to the linearly polarized photon, while p 3 to the circularly polarized photon. The linearly 
polarized beam asymmetry then is given by the cross section differences between polarizing the photons along the Xi 
(<f> = 0°) and 2/j-axis ($ = 90°), which correspond to P 7 along ^fXi, respectively. Similar to Eq. (|28p. we have 

W ± (* = 90°,p)-W l \(* = 0°,p) _ p pjp + 24, 
W ± (<S> = 90°,p) + W ll (<l> = 0°,p) 7 /9 ° + 2p? 1 " 

In terms of the helicity amplitudes the above expression can be written as 

2j - ^ i~ H i-i H u - H l-l H il + H W H iO + H W H 40 
Tjr Tjr rji Tji i Tjr jjr , jji rji 

TJ' r TJ r Tji Tji _1_ TJ' r TJ r ^ TJ^ TJ^ \ 

-•^20-^30 - ^20^30 + H 21 H Z-1 + ^21^3-1/ 

= ^(H\T^ A \H) . (42) 
The above expression is the same as that defined in Ref. 



C. Target polarization asymmetry 

Proceeding to other single polarization observables, we shall derive the angular distributions of the vector meson 
decays into spinless particles and express the polarization observables in terms of the corresponding density matrix 
elements. 
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The polarization status of the target is defined as 

p(N) = ^(lN i +<r-P Ni ) , (43) 

where 1^ is the 2x2 unity matrix in the spin space of the target; Pjv ; is the polarization vector for the target nucleon. 
In the coordinates of the initial states, we define Pjv 4 = -P/Vi (cos <£>, sin$, 0) for polarizing the initial nucleon within 
the Xi-yi plane, and Pjv< = -P/V;(0, 0, A^) for polarizing the initial nucleon along Zi, where is the initial-nucleon 
degree of polarization, and $ is the angle between the polarization vector of the initial nucleon and the production 
plane (a;,-z, plane). 

The polarized target density matrix elements can be related to the elements of the vector meson decay via the 
production amplitudes T: 

A f A-y \i A'. 

The decomposition of the polarization status of the target gives 

p{V) = P° + jZ P kp a , (45) 

a=l 

where p° denotes the density matrix elements with the unpolarized target, and p 1,2 > 3 denotes the elements with the 
target polarized along x%, y% and z^-axis in the initial frame. Explicitly, the elements are given as 



pIk 



A f A-y Ai 



27V 

A f A-y Xi 



,2 



Pkxl 



X f A-y Ai 

E ^ T XvXf,x^Xi T x' v \ f ,x^Xi > ( 46 ) 



2N 

X f A-y Xi 



where A, = ±1 represents the sign of Aj. 

The vector meson decay distribution depends on the properties of the density matrix elements defined above. 
Applying the parity conservation relation and the requirement of the density matrix elements to be Hermitian, we 
obtain the matrix elements for p A A , with a = 1, 2, 3, corresponding to the polarization of the initial nucleon spin 
along Xi, yi, and z%. Substituting the elements into Eq. I|38(l . we obtain the vector meson distribution in terms of the 
different polarization status of the target nucleons: 

W(cos9, </>,$) = W° (cos 9, 0, p° K X[ ) + P m cos &W 1 (cos 6,(1), p{ vK ) 
+P Ni sin 3>W 2 (cos p 2 XvK ) 

+X l P Nt W\cos9,cp,pl vK ) . (47) 



where 



W o (cos9,0,p°) = —[- S m 2 9 + -(3cos 2 9-l)p° 00 



\/2Rep° sin 29 cos cj> - p?_x sin 2 9 cos 2(f)] 



3 

W 1 (cos9, <p, p 1 ) = —[V2Imp{ Q sm29sm(f> + lrnp 1 1 _ 1 sin 2 9sm24>] , 

47T 

W 2 (cos 9, <j),p 2 ) = -^-[p^sin^ + p^cos 2 ^ 

— \/2Repl sin 29 cos <f> — p\_ 1 sin 2 9 cos 20] , 
W 3 (cos 9, (j),p 3 ) = — [V2Rep 3 sin 29 sin 4> + Imp?_ 1 sin 2 9 sin 2S , (48) 

47T 
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where $ again denotes the angle between the polarization vector of the initial nucleon and the production plane. 

As addressed earlier, the corresponding angular distribution of the vector meson decays into pions with the polarized 
target has the same form as Eq. I|H8(I . while the density matrix elements represent different dynamical information. In 
analogy with the polarized photon measurement, the cross sections for polarizing the initial nucleon spin projection 
along ±yi axis, i.e. $ = ±90°, can be obtained by summing all the events together: 

/■7T J i 2lZ P 7T r2lX 

VK T ($ = 90°,p) = / / dnW° (cos 6 7 <p,p°) + P Ni / / dilW 2 (cos 0,(/),p 2 ) 
Je=aJct>=o Je=oJ4>=o 

= W°(p Q )+P Ni W 2 (p 2 ) , (49) 

and 



= -90°, p) = W a (p°) - P Nz W 2 (p 2 ) . (50) 

The polarized target asymmetry is hence defined as the cross section differences between these two polarizations: 

f iy T (<j> = 90°,p)-I^($ = -90°,p) _ p p 2 Q + 2p 2 n 

W t (* = W,p) + W l (* = -W°,p) m P ° m + 2p Q 11 ' 1 ' 

where Pjsr i only depends on the experimental setup. 

Comparing with the polarized beam asymmetry, it shows that the density matrix elements p 1 ' 2,3 defined in the 
target polarization have different properties. In particular, the asymmetry of polarizing the initial nucleon along Xi 
axis would be zero if the vector meson decay events are integrated over 9 and <f>. 

The dynamic significance can be seen more transparently in terms of the helicity amplitudes, i.e. 



T 



TJT TT% TTl TTT \ 



(52) 



which now involves interferences between different helicity amplitudes, and again arrives at the same result as that 
of Ref. @. 



D. Recoil polarization asymmetry 



Similarly, we can also establish the relations for the density matrix elements between the vector meson and the 
recoil nucleon by defining 

p(Nf) = ~(I Nf +a--P Nf ) , (53) 

where is the 2x2 unity matrix in the spin space of the recoil nucleon and Pjy> is the polarization direction 
defined in the final state coordinates, i.e. the frame of (xf, yf, zj) in Fig. ^ In the spin space for the recoil nucleon, 
the decomposition of the recoil polarizations leads to: 

A f X-y Xi 
A f A-y Xi 

PK,\' V = ^/^A 1 ,-A / ,A 1 A I 7 1 ^ A/ A ^ A . 

A f X-y Xi 

PKK = 2iV X! ^/ 7a " a /. a t a 1 t aj,a /: a-,a 1 > (54) 



where A/ = ± 1 represents the sign of A/ . 
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The angular distribution of the vector meson decays can then be expressed as 

W (cos 9, </>,$>) = W° (cos 9, (j), p° Xv X ,J + P Nf cos ^W 1 (cos 9,(1), p\ vK ) 
+P Nf sin $VK 2 (cos 9, cj>, p\ vK ) 

+X f P Nf W 3 (cos 9, 4>, p 3 XvK ) . (55) 

where the W 0,1,2,3 corresponds to distributions with the final state baryon unpolarized, and polarized along Xf, y/, 
and Zf directions. The expressions of W 0, 1,2,3 are the same as Eq. (|48H . 

Similar to the derivation in former subsections, the recoil polarization asymmetry is defined as the cross section 
differences between polarizing the final-state nucleon along ±y/ axis (corresponding to $ = ±90°): 

p _ Wi($ = 90°,p)-W l ($ = -90°,p) _ p Poo + 2p 2 ! ( . 

N f-W r ($ = 90°,p) + W l ($ = -90°,p) Nf p° 00 +2p^ ' l °° J 



which is the familiar form as Eq. (|51() . 

Again, in terms of the helicity amplitudes, the dynamic significance can be seen: 



A 
1 

2 

which is the same as that given by Ref. 



P v i \ ^ r tti Tjr Tjr Tji | rji Tjr jjr jji \ 

N f - Z^^ H 3K H 1X V ~ H 3X V H 1X V + H 4X V H 2X V ~ H iX v H 2xJ 

H^u^H) , (57) 



E. Beam-target double polarization asymmetry 

With the density matrix elements for the photon beams, target nucleon, recoil nucleon, and vector meson, one can 
proceed to the investigation of double polarizations in terms of the vector meson density matrix elements. 

We shall start with the beam-target double polarization, of which experiments are to be available at GRAAL and 
JLab 0. Similar to the single polarizations, the vector meson density matrix elements is related to the polarization 
status of the photon and target nucleon. Namely, we have 

p(V) = Tp( 1 )p(N l )T^ 

= iT(/ 7 + < r-P 7 )(J Wi +a-P Ni )Ti 

= ^T(I 7 I Ni + a- ■ P y I Nt + I J er-P Ni +a- P 7 er • P Ni )T* , (58) 

where one can see that the double polarization only involves the last term in the bracket. The above formula reads 
that the double polarization eventually provides access to the single polarizations as well. 

We introduce the index a and f3 for the polarization status of the photon beams and the target nucleon, and express 
the B-T polarization as 

W (cosB, 4>,p a0 ) = W ao (cos6,(f>,p ao ) 

3 

+ p "W aC ' '(cos 6,4>,p a0 ) 

0=1 

3 

+ ^ P ? P m W af} (cos 0,(t>,p af) ) • (59) 

a, 0=1 

We shall concentrate on the last term, where both the photon and the target nucleon are polarized. The density 
matrix elements thus can be expressed as 

P%K<y) = tf E r A.A / ,X r A 4 pX y X/(7)pf 1 V 1 W)^A / ,^A 5 ■ (60) 
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Given the polarization status of the photon and target nucleon, one can derived the transition elements A , (7) and 
/>[y(iVj) in their spin space. For instance, for the polarization that the photon is linearly polarized along :r-axis 
(a = 1) and the target nucleon is polarized along y-axis (/3 = 2), we have 

Pa^aj, i v ) = X! ^i T ^\ f -\-,-Xi T y v \f,\-,Xi > ( 61 ) 

A / A T A i 

where \\ denotes the sign of A,; . 

The angular distribution again gives access to the experimental measurement of the B-T asymmetry in terms of 
the density matrix elements: 



r jN t _ POO + 2 Pll / fi9 N 
xy n 00 1 o„00 ' V uz 7 



which has the same form as other single polarizations; pjjo and p^i denote the unpolarized density matrix elements 
and Poo + ^Pii = 1 due to the normalization. Undoubtedly, the dynamic information contains in the helicity products 
selected by the polarization status, and again, the B-T polarization experiment will pick up the double polarization 
asymmetry: 

P^Ni * r JT* TT I 7T* TT TT* TT TT* TT 

<~xy ~ 11 3-1 +- H 21" H 4-1 —-#31^1-1 — -^41^2-1 

+H*_ 1 H3i + H^_ x Ha\ — H^^Hn — H^_ 1 H2i 
—Hi H 30 — #20-^40 + H^ H W + Hl H 2 o} 

= -\ E H *axAl<K H bK . (63) 

a,fe, A w , A^ 

which is the same as defined in Ref. Q and and A , are 4x4 and 3x3 Hermitian matrices. 
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FIG. 1: Kinematics for *yNi — » VJV/ in the overall cm. system. 
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FIG. 2: Kinematics for V — > P7 in the overall cm. system where the vector meson has momentum q. The angle 6 C denotes 
the decay direction of the photon, while 6 is the angle between q and the momentum difference between the photon and 
pseudoscalar meson. 



